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Starting from the hydro dynamic equations of binary granular mixtures, we derive an evolution 
equation for the relative velocity of the intruders, which is shown to be coupled to the inertia of 
the smaller particles. The onset of Brazil-nut segregation is explained as a competition between 
the buoyancy and geometric forces: the Archimedean buoyancy force, a buoyancy force due to the 
difference between the energies of two granular species, and two geometric forces, one compressive 
and the other-one tensile in nature, due to the size- difference. We show that inelastic dissipation 
strongly affects the phase diagram of the Brazil nut phenomenon and our model is able to explain 
the experimental results of Breu et al.|l6j|. 

PACS numbers: 45.70.Mg;05.20.Dd 



I. INTRODUCTION 

Segregation is a process in which a homogeneous mix- 
ture of particles of different species becomes spatially 
non-uniform by sorting themselves in terms of their size 
and/or mass 
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Monte Carlo simulations 
of Rosato et al. ;3j clearly demonstrated that the larger 
particles immersed in a sea of smaller particles rise to 
the top when subjected to strong vertical shaking. This 
is the well-known Brazil nut ■phenomenon (BNP). It has 
been explained using the geometrical ideas of percolation 
theory, i.e. in a vibrated-bed the smaller particles are 
more likely to find a void through which they can perco- 
late down to the bottom, leaving the larger-intruder at 
the top0,Q- The arching-eSects Q, whereby the larger 
particle is being supported by the arches of smaller par- 
ticles, can help to assist the percolation-driven segrega- 
tion. The second mechanism of segregation is the convec- 
tive mean flow in the vibrated bed due to the formation 
of convective cells such that the particles move to the 
top through the central-axis |5j- Recently, another mech- 
anism has been proposed, driven by the inertia of the 
intruder [||, which could explain the reverse-buoyancy- 
effect whereupon a light but large particle will sink to 
the bottom of a deep bed under low-frequency shaking. 
In Ref.0, a buoyancy- driven segregation mechanism has 
been proposed, drawing a direct analogy with the stan- 
dard buoyancy forces in a fluid. For other related issues 
on segregation, the reader is referred to the recent review 
article of Rosato et al.0- 

The interplay between size and mass has been con- 
sidered by Hong et aZ.|9| who found that a downward 
intruders' movement occurs as well: reverse Brazil nut 
phenomenon (RBNP). They proposed a phase diagram 
for the BNP/RBNP transition, based on the competition 
between percolation and condensation. Recently, Jenk- 
ins and Yoon investigated the upward downward 



transition employing the hydrodynamic equations for bi- 
nary mixtures. The driving mechanism for segregation 
in the hydrodynamic framework is presumably different 
from that of the percolation-condensation idea and re- 
mains unexplained so far. 

Employing the Enskog-corrected hydrodynamic equa- 
tions for binary mixtures|ll). we investigate the Brazil- 
nut segregation in a dry fluidized granular mixture in 
the absence of bulk convection. The purpose of this pa- 
per is three-fold: firstly, to derive a time-evolution equa- 
tion for the relative velocity of a single intruder, taking 
into account the non-equipartition of granular energy; 
secondly, to explain the driving mechanism for Brazil- 
nut segregation in terms of the buoyancy and geometric 
forces. Lastly, based on a simple model for energy non- 
equipartition, we will show how the inelastic dissipation 
determines the regimes of BNP and RBNP. 



II. HYDRODYNAMICS OF GRANULAR 
MIXTURES 

The validity of the hydrodynamic approach even in 
the dense granular flows has recently been justified via 
the comparison of theory with various experiments|l2|- 
here one has to be careful in choosing the appropri- 
ate constitutive model for pressure, viscosity, dissipa- 
tion, etc. The constitutive model that we have used has 
been validated by performing MD simulations of binary 
mixtures [l3|. We consider a binary mixture of slightly 
inelastic, smooth particles (disks/spheres) with radii 
(i = I, s, where index I stands for large and s for small), 
mass mi and number density m. The species mass den- 
sity is @j(x,i) = miUi = Pi4>i, where pi is the material 
density of species i and <pi is its volume fraction. The 
total mass-density, g(x, t), and the total number-density, 
n(x, t), are just the sums over their respective species val- 
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ues. The dissipative nature of particle collisions is taken 
into account through the normal coefficient of restitution 
e.y, with e,j = eji and < < 1. 

Assuming unidirectional flow (iij = (0, Vi(y, t), 0); 
d/dx = 0, d/dz = and <9/ch/ 7^ 0) and neglecting vis- 
cous stresses, the momentum balance equation for species 
i can be written as 



dvj 
1 dt 



dpi 
dy 



(1) 



Pi 



UiZiTi, with Zj = 1 + JQ. 



(2) 



Here is the compressibility factor of species i and 
A'y = (^jS'y (1 + Rij) d /2, with being the radial distri- 
bution function at contact and Rij = Tijtj the size-ratio. 
Note that the ify are related to the collisional compo- 
nent of the partial pressure, having a weak dependence 
on inelasticity which we neglect, and Zj — > 1 in the dilute 
limit </> — > 0. We shall return back to discuss more about 
the compressibility factor later. 

After some algebraic manipulations with the momen- 
tum balance equations and the equation of state, we ob- 
tain the following evolution equation for the relative ve- 
locity of the larger particles, v\ = vi — v s , 
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An explicit expression for the momentum source term, 
Ti , can be obtained using the Maxwellian velocity distri- 
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(Note that the mass balance equations are identically sat- 
isfied for unidirectional flows.) Here pi is the partial pres- 
sure of species i, and g the gravitational acceleration act- 
ing along the negative y-axis; Tj is the momentum source 
term which arises solely due to the interactions between 
unlike particles and YLi=i a Ft = [ijj. The assumption mi ~M 
of negligible viscous stresses is justified if there is no over- 
all mean flow in the system, or if the spatial variation of 
Vi(y, t) is small. 

To obtain constitutive relations for partial pressures 
we take into account the breakdown of equipartition 
of energy between the two species (in the equation of 
state) as found in many recent theoretical and numerical 
studiesl7l ll ll Il3t Il4j and also confirmed in vibrofluidized 
experiments [T3 • We assume that the single particle ve- 
locity distribution function of species i is a Maxwellian at 
its own granular energy Ti, where Tj = ^ < Cj-Cj >, 
with d = 2 and 3 for disks and spheres, respectively, 
Cj = Cj — u being the peculiar velocity, Cj the instan- 
taneous particle velocity and u = g^ 1 J2i=i s ft u » the 
mixture velocity. The equation of state for the partial 
pressure of species i can then be written as: 



where T = n _1 J2 i=l s niTi = J2 i=l s ^jTj is the mixture 
granular energy, £j — rii/n the number-fraction of species 
i, mis = m; + m s and ri s = r\ + r s . With additional as- 
sumptions of weak gradients in species number densities 
and granular energy, and retaining terms of the same or- 
der in the single intruder limit (ni << n s ), the evolution 
equation can be considerably simplified to 
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This is our time-evolution equation for the relative veloc- 
ity of a single intruder: the first term on the right hand 
side is the net gravitational force acting on the intruder, 
the second term is a 'Stokesian-like' drag force and the 
third term represents a weighted coupling with the inertia 
of the smaller particles. It is interesting to recall the work 
of Shinbrot and Muzziojg who argued that the onset of 
reverse-buoyancy would crucially depend on the inertia 
of the smaller particles- a detailed analysis of eq. (0 with 
appropriate boundary conditions is left out for a future 
investigation. In this paper, we are only interested in the 
steady-state solution of the above equation. In typical 
situations where one can neglect the last term (e.g. if 
the intruder is much heavier than the smaller particles), 
we end up with the familiar evolution equation where 
the inertia of the intruder is being balanced by the net 
gravitational force and the drag force. Only in this case, 
the interplay between the gravitational and drag force 
will eventually decide whether the intruder rises or sinks. 
Neglecting transient effects, the steady relative velocity 
of the intruder can be obtained from 
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Setting this relative velocity to zero, we obtain the crite- 
rion for the transition from BNP to RBNP: 



ZsT s 



mi = 0, 



(J) 



which agrees with the expression of Jenkins & Yoon[T(| 
for the case of equal granular energies (T; — T s ). As 
such, it is not evident from this expression what the driv- 
ing mechanism for segregation is. Thus, we need to an- 
swer several questions. Can we recast the segregation 
criterion in terms of the well-known Archimedean and 
thermal buoyancy forces? Is there any new force, and 
what could be the physical origin of such forces? 
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III. DRIVING MECHANISM: SEGREGATION 
FORCES 

To understand the origin of segregation in the present 
framework, we now decompose the net gravitational force 
in eq. (J5J) for a single intruder in the following manner: 
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where Vi is the volume of a particle of species i. The first 
term, Fg = Vi(p s — p/)g, is the effective Archimedean 
buoyancy force which arises due to the weight of the dis- 
placed volume of the intruder (VI). The second term, 
Fg oc (T; — T s ), represents the buoyancy force due to 
the difference between the two species granular energies. 
This, being an analog of the thermal buoyancy, may be 
termed the pseudo-thermal buoyancy force. 

There are two more terms in eq. JHJ which do not ap- 
pear to be related to standard buoyancy arguments. The 
third term is negative definite, and vanishes identically if 
the intruder and the smaller particles have the same size. 
Note that ej* = (Vi/V s — 1) is the volumetric strain. Thus, 
Fg* = —m s ge s * is a static compressive force to overcome 
the barrier of the compressive volumetric strain arising 
out of the size- disparity between the intruder and the 
smaller particles. 

The fourth term in eq. JSJ, oc (Zi/Z s — 1), van- 
ishes in the dilute limit (f> — > 0. It can be verified 
that (ZijZ s — 1) also vanishes identically, irrespective 
of the total volume fraction, if the particles are of the 
same size. (Note that we have neglected the weak- 
dependence of Zi on inelasticity.) Thus, the origin of 
this force is also tied to the size- disparity as in the third 
term F g ^. An interesting physical interpretation can be 
made if we consider the dense limit with a single in- 
truder (0, « fa): (Zi/Z s - 1) oc Rf s for R ls » 1. 
Hence e^ y ™ = (Zi/Z s — 1) > can be associated with 
a weighted volumetric strain, tensile in nature. Thus, 
Fge n = m sE e i vn is a dynamic tensile force that arises 
from the excess pressure difference due to the nonideal 
(collisional) interactions between the intruder and the 
displaced smaller particles. 

Thus, the geometric effects due to the size- disparity 
contribute two new types of segregation forces: 
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the former is a static, compressive force and the latter 
is a dynamic, tensile force. On the whole, the collisional 
interactions help to reduce the net compressive force that 
the intruder has to overcome. 

A question naturally arises as to whether we could 
get back the standard Archimedes law from eq. (JSJ if 
we take the corresponding fluid limit, i.e. a large par- 
ticle being immersed in a sea of small particles with 
n >> r s . In this limit it immediately follows that 



m s (Vi/V s -l) = -Ff e and hence F ge = 0. Thus, 
the net gravitational force on a particle falling/rising in 
an otherwise quiscent fluid (at the same temperature) is 
nothing but the standard Archimedean buoyancy force, 
F = Fg = g(p s — pi)Vi. It is worth recalling that when 
there is no size-disparity (n = r s ), the geometric forces 
are identically zero. Hence the behaviour of a heavier 
particle in a sea of equal-size lighter particles is similar 
to that of a particle in a fluid. 

To clarify our segregation mechanism, we show the 
variations of different segregation forces with the size- 
ratio in Fig. n f° r t ne two-dimensional case of equal 
density particles (pi — p s ) in the single intruder limit 
(4>i/4>s = 10~ 8 ) at a total solid fraction of <p = 0.7, with 
the restitution coefficient being set to 0.9. For illustra- 
tive purposes, we have calculated the energy ratio, Ti/T s , 
(see the lower inset in Fig. ^) from the model of Barrat 
and Trizac|l4j. For this case, the Archimedean buoyancy 
force is identically zero, and the total geometric force re- 
mains negative, as seen from the upper inset in Fig. ^ 
The pseudo-thermal buoyancy force is, however, positive. 
Thus, the competition between the pseudo-thermal buoy- 
ancy force and the geometric forces leads to a critical size- 
ratio above which the intruder will rise for this case. (For 
the corresponding purely elastic case (e = 1 and Fg = 0) , 
the net force is F = F ge < and hence the larger particle 
will sink to the bottom.) This mechanism holds also for 
the more general case (pi ^ p s and Fg ^ 0) for which 
the total buoyancy force (Fb = Fg +Fg) competes with 
the geometric forces (F ge = Fj£ + Fg^ n ) to determine 
the transition from BNP to RBNP; the inclusion of dissi- 
pation merely affects the location of the transition point 
(see Fig. |21 and the discussion below, for details). 



IV. PHASE DIAGRAM AND DISCUSSION 

A typical phase diagram in the single intruder limit 
(4>i/(j> s = 10 -8 ), delineating the regimes between BNP 
and RBNP, is shown in Fig. [21 for the two-dimensional 
case, with other parameters as in Fig. ^ (The quali- 
tative features of the corresponding phase diagram for 
the three-dimensional case are similar.) Focussing on 
the purely elastic case (e = 1), we note that a transi- 
tion from BNP to RBNP can occur following two paths 
(denoted by two arrows), one along the constant mass- 
ratio with decreasing size-ratio and the other along the 
constant size-ratio with increasing mass ratio. In both 
cases, the Archimedean buoyancy force balances the net 
geometric forces at the transition point. 

Comparison between the elastic (e = 1) and inelas- 
tic (e = 0.9) cases in Fig. [21 clearly shows that the 
non-equipartition of granular energy, responsible for the 
pseudo-thermal buoyancy force Fg, has a dramatic ef- 
fect in reducing the regime of RBNP, and decreasing the 
value of e reduces the size of this regime further. For the 
case of a mixture with equal volume fractions (<f>i — <f> s ), 
however, the regime of RBNP is much larger as seen from 
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FIG. 1: Variations of segregation forces (F/m s g) with the 
size-ratio for pi = p s at e = 0.9; see text for other details. 
The upper inset shows the corresponding static and dynamic 
contributions to the total geometric force (F ge = Fgl + Fg^ n ). 
The lower inset shows the variation of Ti/T a with the size- 
ratio 0. 




FIG. 2: Phase diagram for BNP/RBNP in two dimensions: 
4> = 0.7 and (f>i/4>s = 10 -8 . Left inset: phase diagram with 
e = 0.9, 4>i/4> s = 10 -8 (solid curve) and <f>i/<j> s = 1 (dashed 
curve). Right inset: phase diagram with e = 0.9, (j>i/<f>s = 1, 
4> — 0.7 (dashed curve) and (j> — 0.4 (dot-dashed curve). 

the left-inset of Fig. [21 This observation is in qualitative 
agreement with the recent experimental results of Brcu 
et q/.|l6j| who found that the 1 reverse Brazil-nut effect is 
completely destroyed if <j) a » 4>i ■ 

The right-inset of Fig. [21 shows that the size of 



the RBNP-regime also increases with decreasing over- 
all mean volume fraction. Since in vibrated-bed experi- 
ments increasing the shaking amplitude is equivalent to 
decreasing the mean volume fraction, our observation ex- 
plains another interesting result of Breu et aL[l6( that for 
a given mixture with specified size- and mass-ratio, the 
final state is that of RBNP at sufficiently high accelera- 
tions (see Fig. [21 in Ref. J16j|). 

We need to point out that calculating the energy ratio 
(Ti/T s , we made the assumption that e%j — e. Using 
a variable restitution coefficient, our phase-diagram at 
large mass-ratios will be modified, but the proposed seg- 
regation mechanism and the qualitative features of the 
phase-diagram remain intact. Even though the model of 
Barrat & Trizac^4| is strictly valid for a homogeneous 
mixture with stochastic-driving, it has recently been 
verified in vibrofluidized-bed experiments under strong 
shaking 

To compare our segregation mechanism with others, 
we note that the scaling of the geometric forces (ex Rf s ) 
suggests that they can be compared to the effective per- 
colation force of Rosato et al.y|, and hence we have a 
competition between buoyancy and percolation forces. 
In the percolation-condensation mechanism of Hong et 
al.0, the condensation is driven by the two-species hav- 
ing different energies. If one equates their driving force 
due to condensation-tendency with an effective buoyancy 
force, then our mechanism could be equivalent to that of 
Hong et al. However, there is no direct one-to-one anal- 
ogy between our hydrodynamic segregation mechanism 
and the percolation-condensation mechanism. 

In conclusion, we have identified four different types 
of segregation forces: apart from the Archimedean buoy- 
ancy force and an analog of the thermal buoyancy force, 
there are two additional forces, the origin of both is tied 
to the size-disparity between the intruder and the smaller 
particles. We have demonstrated that the competition 
between the buoyancy and geometric forces determines 
the onset of segregation in the present scenario, and the 
inclusion of the pseudo-thermal buoyancy force (due to 
inelastic dissipation) further enhances the possibility of 
BNP. While the possibility of RBNP is rather limited in 
the single intruder limit, even at moderate dissipation- 
levels, either increasing the relative volume fraction of 
the intruders or decreasing the mean volume fraction 
enhances its likeliness as in the experiments of Breu et 
al.p3. 
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